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Abstract 

Mh 

^ Non-radiative coupling between conductive coils is a candidate mechanism for wireless energy 

transfer applications. In this paper, we propose a power relay system based on a near-field metama- 
p ^ terial superlens, and present a thorough theoretical analysis of this system. We use time-harmonic 

^ circuit formalism to describe all interactions between two coils attached to external circuits and 

^ coils are found in the point-dipole approximation using Sommerfeld integrals, which are reduced 

^ to standard special functions in the long- wavelength limit. We show that, even with a realistic 

magnetic loss tangent of order 0.1, the power transfer efficiency with the slab can be an order of 
magnitude greater than free-space efficiency when the load resistance exceeds a certain threshold 

(N 

K»" value. We also find that the volume occupied by the metamaterial between the coils can be greatly 

^ compressed by employing magnetic permeability with a large anisotropy ratio. 

(N 

<N 
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I. INTRODUCTION 



The explosive growth in the use of cordless hand-held electronic devices, electrical vehicles 
and appliances in the recent years is stimulating interest in wireless power source^^. Wire- 
less communications between stationary and mobile devices make it possible to remotely 
control devices, giving them a certain degree of autonomy. However, since all electrical de- 
vices depend on an energy supply, it must be either carried with the device, or harvested 
from the environment as needed. 

Energy harvesting can include mechanical, thermal, chemical, gravitational or electro- 
magnetic energy, although only the latter is suitable for most applications. Radiative trans- 
fer of electromagnetic energy is limited by absorption and scattering in the atmosphere and 
requires a direct line of sight between the source and the device^. At high intensities, it 
also presents a challenging electromagnetic interference problem. While in the radiation 
flux the magnitudes of electric and magnetic fields are always comparable, in the near field 
of magnetic dipole antennas the EjH ratio can be strongly suppressed, thus reducing the 
interaction with biological and other environmental objects whose properties are almost 
purely dielectric. This bio-friendly behavior of inductively coupled resonators is widely used 
in RFID transponder^. Thus, non-radiative inductive coupling between high-frequency 
circuits within their near-field zone is an attractive option for wireless power transfer appli- 
cations. 

Although inductive near-field coupling between small coils of radius r ^ A separated 
by a distance d such that r ^ (i ^ A is generally small - proportional to {r / X){r / dY - 
it can be resonantly enhanced if the coils are either self-resonant or connected to external 
resonant circuits. Wireless power transfer over a distance of 8 times the radius of the coils 
with efficiency 40% has been experimentally demonstrated by Soljacic et al.'^. However, the 
power transfer efficiency of their device drops steeply as a function of the distance between 
the coils, as well as a function of the resistive load attached to either of the coils. To address 
these issues, we propose a relay system based on the concept of the near-field superlend^, 
which can greatly enhance both the range of distances and the load resistance at which the 
power transfer efficiency is high enough to be practically useful. 

Shortly after the discovery of the negative-index superlen^^^II]^ \i realized that its 
superresolution is related to the enhancement of the transfer function of the near fields, in 
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addition to its aberration-free refractive properties^^"^. Reduced near-field-only superlenses 
based on negative permittivity media were theoretically proposecP^HSl implemented 
shortly afterward j^°*^^l An important aspect of negative-e lenses is that they enhance trans- 
mission of only TM-polarized waves (with magnetic field normal to the plane of propagation). 
In two-dimensional geometries, TM waves can be generated by either in-plane polarized elec- 
tric dipoles or out-of-plane magnetic line currents. The near-perfect transfer of magnetic, 
as well as electric, fields has been demonstrated in two-dimensional systems with purely di- 
electric superlense4^^^. Similar results are expected with purely magnetic superlense^^^^ 
due to electromagnetic duality. 

In contrast with their two-dimensional counterparts, three-dimensional magnetic point 
dipoles are difficult to couple using a dielectric superlens, because the TM wave component 
of their magnetic field is suppressed in the near-field relative to the TE wave component by a 
small factor of order (rf/A)^ (see Refs.'^^'^ and equation 24 below). Fortunately, for applica- 



tions relying on radio or microwave frequencies, artificial media with negative magnetic per- 
meability have been designed and demonstratecP^'^^. Negative-// metamaterial lenses have 
been proposed for improved magnetic resonance imagin^^^M) g^^^^ wireless energy transfei'^^l 
applications. In this article, we focus on finding the optimum coupling regime that maxi- 
mizes the amount of energy transferred to the useful load. In performing the analysis, we 
allow the metamaterial to be anisotropic. Earlier studies suggested that strongly anisotropic 
media can serve as enhanced imaging systems and may dramatically enhance the transfer 
function!I2I2BH3Q] Eliminating the need for strong magnetic response in one or more directions 
is potentially interesting as it could reduce the complexity and cost of such metamaterials. 
However, our findings indicate that the optimum coupling is obtained in the regime where all 
three components of fi tensor are simultaneously negative. On the other hand, simultaneous 
control over all three principal values of /i enables reduction of the superlens thicknes^^^'^. 

This paper is organized as follows. First, we introduce a circuit model of the induc- 
tive coupling between two coils and one metamaterial slab, and define the power transfer 
efficiency equivalent to the one introduced in Ref.'^. Then, we calculate the mutual and 
self-inductances of the coils in the presence of the slab, and show that a lossy metamate- 
rial introduces resistive contributions to the otherwise purely reactive response of perfectly 
conducting coils. The coils are described as point dipoles where applicable, and the metama- 
terial is replaced with an equivalent homogeneous, possibly anisotropic, medium. Then, we 




(a) (b) 

FIG. 1. (Color online) Left: The schematic of the device consisting of two single-turn coils coupled 
through a slab of homogeneous, possibly anisotropic medium, of thickness D. The z axis is normal 
to the slab, with the origin (0) on the left face of the slab. The coils are positioned on the z axis, 
with their magnetic dipole moments oriented in the x-direction. The distance from the center of 
coil 1 to the left face of the slab is di, and from coil 2 to the right face of the slab it is d2- Right: 
Circuit diagram of the device. 

analyze the power transfer efficiency and find the optimum couphng regime. Finally, several 
important aspects of this relay system are discussed: the effect of the dipole orientation 
relative to the slab, contribution of the TM waves and the possibility of enhancing magnetic 
fields with a dielectric-only metamaterial, and the role of anisotropy. 



II. THE CIRCUIT MODEL OF MAGNETIC DIPOLE COUPLING 

Near-field coupling between magnetic dipoles is an attractive option for radiation-less 
power transfer. In contrast with electric dipoles, magnetic field in the vicinity of magnetic 
dipoles dominates strongly over the electric field; thus, magnetic dipole-based energy transfer 
system is expected to interfere less with the function of other devices and live biological 
objects, including humans. A simple implementation of a magnetic dipole - a circular, 
planar, single-turn loop of wire - is assumed in this study (see Fig. [l](a)). The magnetic 
moment of a wire loop relates to its current i and area A = ttB? as mo = lAn, where fi is a 
unit vector normal to the plane of the loop. 

Electromagnetic interaction between two small wire loops can be described with the aid 
of the coupled-coil theorj^, which states that the magnetic fiux through coil m is expressible 
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through the currents in all participating coils: 



^ (-SjTi ■ TL^j^ Ayji ^ ^ LYYin^ni (l) 

71=1 

where is the area of the m-th coil, Bm is magnetic flux density at the center of the loop, 
and the coefficients Lmn are called self-inductances for m = n, and mutual inductances for 
m n. The radius of each coil is assumed to be much smaller than both the free-space 
wavelength and the spatial scale of variation of the magnetic fields generated by all other 
coils. The expressions ([T]) remain valid even if the coils are located in the vicinity of objects 
with strong electromagnetic response, for example, a metamaterial slab, as shown in Fig. [1} 
Inductances Lmn = Lnm are affected by the coils environment; they are calculated for the 
two-coil-and-one-slab system below. 

Using Faraday's law of induction, the relationships ([T]) can be re-written in terms of the 
induced electromotive force (EMF), i.e., inductively induced voltage in each coil: 

di ^ 

'^m ~ ^ ^ ~ ^ ^ Lmn'i'n- (2) 

n=l n=l 

We assume time-harmonic dependence of the currents and fields as 2 oc e~^'^^. Assume that 
coils 1 and 2 are connected to external circuits with resistances Rm and capacitances Cm, 
m = 1,2; additionally, the circuit of coil 1 includes a voltage generator with a fixed external 
EMF S. From the two voltage balance equations, 

Vi = £ + iiRi + = -juj{Liiii + Li2i2), 

V2 = + —-^ = -j(^iL22i2 + L2lil), (3) 
JUU2 

the currents in the two coils are expressed through the external EMF, giving 

-S 



eff ' 



■ L21 . . L21S , , 

^2 = -juj—ti = jw-^g— , (4) 

Zj2 Zj^ ^2 



where we have introduced effective coil impedances 



^1 — ^1 + a; — - — , 
^2 

Zm = Rm^3^{Lmm ^-^),m=l,2. (5) 

IjJ Urn 
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The power dissipated in each circuit is given by the sum of terms Pm — |Re(V^i^) = 
|Re(Z^)|i^P, summed over all circuit elements with non-vanishing resistance; explicitly, 
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where we have introduced effective resistances 

Rf = ReZf = i?i - wImLii + uHm^^^^, (7) 

Rf = ReZ2 = i?2 - a;ImL22 + w^Im^^^^ . (8) 

Note that the self-inductanccs of coils in the presence of a lossy slab may have a non-zero 
imaginary part, which accounts for the resistive, dielectric and magnetic losses in the slab. 
In the following, we will be neglecting the terms quadratic in frequency, which originate 
from the resistive part of mutual impedance Z12 = jujLi2. Because of this resistive part, 
one should distinguish the power dissipated in the useful resistive load R2, 

P'2 = \R2\i2\\ (9) 

from the total power dissipated in circuit 2 (^2)- 

The fraction of power delivered to the resistive load R2 relative to the total consumed 
power is 

P2 R2 X 

where 

p peff, ,21 r 12 

P2^ R2^ IL21I 

In order to make the power transfer efficiency r] close to 100%, one wants x to be as large 
as possible, and the ratio Rf /R2 > 1 as small as possible. To maximize x, it is beneficial 
to minimize \Z2\ by eliminating ImZ2; this is achieved by tuning circuit 2 to its resonance 
by adjusting the capacitance C2, such that ReL22 — ^^^^^ — 0. In the following, we assume 
that Z2 — Rf, and thus, 

/,2|7- 12 

_ ^ 1-^21 1 f^r,\ 

^ ~ RfRf ■ ^ ^ 

We observe that the power transfer efficiency rj is not sensitive to the magnitude of Zf, 
but only to its resistive part. That is because the ratio x, as well as efficiency rj, do not 
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reflect the actual power extracted from the voltage generator £ by the load. In order to 
characterize the effect of quenching of the source in the presence of the slab, it is convenient 
to introduce an additional figure of merit, Xid) as the ratio of the useful transmitted power 
to the power extracted by an "ideal" load connected directly to the voltage source, i.e. 

Pi, 





2 


L2I 


2 p2 




ZfZ2 


2 



f R2U; 


-^21 




Uf 1 


zf\ 





_ ^ 1-^211 -'^■2 / -n.2 ^\^2l\ \ 

Xid — I ryoff rv 1 9 ^ 1 nr-ff I ^cff I ' V-'^l 



2R2 ■ 



where in the latter limit we again assumed \Z2\ ^ Rf . Note that both x and Xtd may be 
greater than one, whereas r] <lhy definition. 

In order to maximize Xid, the source coil must operate at a frequency at which \Z^\ is 
minimal. By varying the capacitance Ci, one can eliminate the imaginary (reactive) part of 



Zl , in that regime, the figure of merit (13) becomes 
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_ / 0^\L2l\R2 \ ,x 
Xid — \ peff peff J " \^^) 



We have expressed two useful figures of merit through the modified coil impedances, 
Zmn = jojLmn, which are calculated in the subsequent sections. It is interesting to see 
how they relate to the shift in resonance frequency induced by the slab. The resonance 
frequencies are obtained by equating Zf = 0, which gives a polynomial equation for the 
frequency: 

ZiZ2 + u^LuL2i=0. (15) 

Neglecting resistive loads (i?i = i?2 = 0), a biquadratic equation is obtained for the reso- 
nance frequency: 

co^ det{L„,n)CiC2 - uj\L^iC^ + L22C2) + 1 = 0. (16) 

Its solutions are 



'^^ - 2{ujfuf - Q-^) ' ^^^^ 

where cu^ = 1 / {CmLmm) and f2~^ = Li2L2iCiC2- For the frequency difference we obtain 

For a symmetric system (wi = U2 = ojq), the frequency shift to the lowest order in is 
Au = u)+ — ~ Q^'^Uq. Since oc L12-L21 oc x Xid, we can see that in the small 
coupling {x 1)) lossless limit, a linear relationship exists between the power coupling 



,2 ,a, v^^. ■ (18) 



efficiencies and (Aw)^. This relationship is however hnear only if the slab contributions to 
self-inductances 22, and consequently to the resonance frequencies Ui^2, are disregarded. 

III. CALCULATION OF MUTUAL INDUCTANCE IN THE POINT-DIPOLE AP- 
PROXIMATION 

To describe the operation of current- carrying coils in the presence of a metamaterial slab, 
one needs to know the inductances Lmn- Self-inductance of a wire loop is a function of the 
loop radius as well as the wire thickness. On the other hand, mutual inductance L12 can 
be determined with good accuracy from the value of magnetic flux density B generated by 
coil 1 evaluated at the center of coil 2. The distribution of B in space can be found in the 
point-dipole approximation, i.e. by neglecting the size of the source coil. In this section, we 
calculate L12 assuming that both coils are negligibly small in comparison with the spatial 
scale of the magnetic field gradient. The problem is then reduced to the well-known problem 
of a dipole antenna over a fiat layer of homogeneous medium. For an infinitely extended 
substrate, this problem was solved by SommerfelcP^, and the solution is often referred to 
as the Sommerfeld integral. For an infinitely wide slab of finite thickness, the solutions 
in terms of similar integrals were reported by Weng CheW^. Here, we generalize these 
solutions for a slab with uniaxially anisotropic properties. In the quasistatic limit, we are 
able to obtain closed-form expressions for all fields in terms of a well-known special function, 
Lerch transcendent 

We assume that a slab of thickness D, which occupies the space between < z < D, is 
embedded into homogeneous medium, such as air, with relative permittivity e„ and relative 
permeability fi^. The relative permittivity and permeability tensors of the slab are assumed 
uniaxial, with Cartesian components 

e = diag(e^, ey = e^;, e^),^ = diag(yUa;, ^-y = ^x, fJ'z)- (19) 

Since the principal axis of the material properties is aligned with the normal to the slab, 
the cylindrical symmetry of the problem is preserved, and one can use the approach of 
even with the medium anisotropy; the benefits of anisotropy are explained below. 
The source, a point magnetic dipole with magnetic moment mi = iiAi, oriented parallel to 
the slab, is placed on the z-axis at z = —di. 
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On the reflection side of tlie slab, in the region with — lii < 2; < 0, the axial components 
of electric and magnetic field can be written ad^ 



mi 



cos I 



WyU^,/ioml 

87 



/oo 
-00 

/oo f^2 
dkpfHi'\k,p)e^''^^' [e^"^' + RTM{kp)e-^^^'] 
00 '^z 



(20) 



where kz = yk^eyfi^ — k'^, and Rte and Rtm are reflection coefficients for the electric field 
in TE and TM polarizations, which are derived below. On the transmission side {z > D), 
the fields are 



mi 



E. 



Stt 
up^Pomi 



.2f7(l)n. r.\rr^_(U \Jk,(z-D) 



COS0 / dkpk Hl '{kpp)TTE{kp)e^ 



sm I 



67r 



k 



dkp-^H[^\kpp)TTM{kp)e 



jk4z-D) 



(21) 



To facilitate the derivations, introduce the spectral components of all fields according to 
RefP 



H{p,^,z) 
E{p, 0, z) 



dkpE{kp, z, 



(22) 



and notice that the problem reduces to the one-dimensional propagation of plane waves 
E[kp, z, 0) in the ^-direction through the layered medium. The reflection and transmission 



coefficients, calculated in accordance with our phase normalization choices seen from (21), 
are calculated in the Appendix; see equations ( 72|[77 ). 

The remaining four components of the E and H vector fields, are expressed through 
Hz{kp,z) and Ez{kp,z) using the formulas of Chew-^ (in Cartesian coordinates): 



Hs{kp,(j),z) 



Vs— \- joje^eoz X VsE^ 



Es{kp,(l),z) 



1 

kl 



dz 

dE^ 
dz 



jup^poz X VsHz 



(23) 



where subscript s indicates the transverse part of a vector. In particular, for Hj. values on 
the axis, in the region z > D (where the second coil would be placed), we obtain 



H^{kp,p = 0, 



jmi 
IGvr 



kpkzTTE 



Jk,{z-D) 



(24) 
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where k'^ = e^fivk^. In deriving equation (24), we have used J[{0) = 1/2 and the oddness of 



the Bessel function Yi{kpp) = lmH[^\kpp). 

So far, the calculations have been exact. To simplify further analysis, assume that the 
distance between the coils is deeply sub-wavelength, i.e. d = D + di + d2 Xv = 2Ti/ky. 
Since we are interested only in the near- field region around it, only Fourier components with 
\kp\ 3> /c„ are important, and the exact formulas above can be greatly simplified. Under this 
assumption, the following approximations can be made (as in Ref.^^J): 

Xo = -jK ~ \kp\, (25) 



Xte = \J k^pfJ'x/fJ'z - l^'x^xkl ^ l^pl V/I^7a^5 (26) 
Xtm = J kje^/e^ - e^fi^k'^ ^ \kp\^/eje~,. (27) 



Note that we implicitly assume that signRe/Xa; = signRe/i^ and signRee^; = signRee^; in 
these regimes, spectral components evanescent in vacuum regions are also evanescent (or 
exponentially growing) inside the slab. The reflection and transmission functions in this 
quasistatic approximation are 



-TE 



aj,^e°'TE\kp\D _^ lj^^(,-aTE\kp\D- 



(28) 



where arE ^ V^W^, arE ^ - + i) ^ ^te ^ - j-) , Qte = - and 

hxE = are material constants; analogous expressions for TM components are obtained 

by the substitution fi ^ e. 

Magnetic field intensity at the position of the second dipole becomes 

H^ip = 0,z = D + d2) = ^ dkp [klTTEi\kp\) + klTTMi\kp\)] e'l'^"!'^^ 

POO 

/ dkp[klTTEikp) + klTTMikp)]e-''''^\ (29) 
Jo 



Stt 

Similarly, at the position of the first dipole we have 



H,{p = 0,z = -di) = ^ rdkpkl [{1 + RTEikp)e-^''"'^) + kl{l + RTMikp)e-"'^''^J^O) 



The integral in (30) diverges, as it includes the field generated by a point dipole. The 
reflected portion of this field is nevertheless finite, and it will be used in the next section for 
a calculation of slab effect on self-impedance. 
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We observe that the contribution of the TM wave is suppressed by a very small factor, of 
order (koD)^. Practically, it means that dielectric slabs cannot provide strong enhancement 
of the mutual inductance, thus the need for magnetic metamaterial slabs. 



The only remaining integral in equation (29) can be reduced to the standard Lerch 
transcendent function using the substitution t = exp{—2aDkp). Noticing that 

'b, dt = 2$i(— ,s + (31) 

a 

we obtain the flux through coil 2: 

$2 = B^ip = 0,z = D + d,)A2 = -'^'^'f^' r dk, [klTTE + klTTM] e-'^'^ 

8vr Jo 

—^-v^-oilMM \hTE 1 ^ / bTE o N I ^TM kl BtM . Ar,\ 
= Z 77, ,3,Mte) + J7. ,1,mtm13?) 

where ute,tm = {o:te,tmD + di + d2)/{2aTE,TMD) > 0. In the following, we will be 
neglecting the TM contribution, and omitting the TE subscript for all variables. The mutual 



inductance coefficient is obtained from expression (32) by dividing out the current i^. 

= = ^^(2^^^(-a'''")- 

For positive real u, the function ^l{'j,3,u) is analytic in both 7 and u, except the line 
7 = [1; +00) where it has a branch cut discontinuity in 7. Since the branch cut discontinuity 
affects only the phase of and the physical result depends only on we can treat this 
function as continuous. 

To highlight the effect of the slab on the inductive coupling, one should compare L21 



given by (33) to its value in the absence of the slab, i.e. when /ij, = /i^ = 1 and fiy = 1. 



Substituting a = h = —4 and b = g = into (33) and using $l(0, 3, m) = m , we obtain 

^21 ^ (^<^>L(0,3,rf/(2D)) = 4^^3 • (34) 

where d = D + di + d2 is the distance between the coils. 



IV. THE EFFECT OF THE SLAB ON SELF-INDUCTANCE 



The contribution of the slab response to the self-inductance of coil 1 can be also estimated 
in the point-dipole approximation using the following modeP^. Consider a point magnetic 
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dipole fhi = rriix excited by the sum of external magnetic field, , and the field reflected 
from the slab, 

= G7jm,, (35) 



where 



1 /"OO 

firef _ 

XX Q„ 



(36) 



is the reflected portion of the Green's function, as seen from equation (30). The total 
magnetic moment equals 

m, = aUHr + H:^f), (37) 

where is the intrinsic magnetic polarizability of the dipole unperturbed by the slab. 
Resolving equations ( 35p7 ), we find rrii = 0^/(1 — o.mG'^xx)-^T^ ^ which means that the 
effective polarizability of the dipole in the presence of the slab is 



[a. 



-1 _ ^ref\-l 
a ^xx J 



(38) 



From the relationship $i = Luii = Luirii/Ai we find that this polarizability relates to 
self-inductance as Ln = fiyfioAl/am', thus. 



(39) 



where L^i^ is the self-inductance of coil 1 without the slab. In the quasistatic approximation, 
and neglecting the weak contribution of the TM wave, we obtain the slab contribution to 
self-inductance Ln: 
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r{l) _ _,. _ f^vf^oA 9 



^l{--,S,v)-^l{--,3,v + 1) 
a a 



(40) 



where v = di/ (aD). Using the identity 



1 



$l(7, s,u + 1) = - ($^(7, s,u)-u ') 



(41) 



expression (40) for f > can be reduced to 

(1) _ /it,/io^? g 



'11 



47r(2aD)3a 



a\ ^ , b . a 
+ t)'^l{—,S,v)--v- 
0/ a 



(42) 



The expression for L22 is obtained by replacing v = di/{aD) in (42) with w = d2/{aD) 
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V. POWER TRANSFER EFFICIENCY ANALYSIS 



Now we turn to the most reahstic scenario where the self-resonant coils are situated in 
free space (e„ = fiy = 1). Optimization of the power transfer efficiency, 

R,-ulm[l<^L{-b/a,3,w)]l + x' ^ ' 

where 

\a) (fli-a)Im[;4i(-''/a,3,t')])(fl2-'^Im[;*i(-''/a,3,«')l) 
requires finding the optimum balance between the enhancement of L12 and the growth of 

ImL^^''22; it is easy to see that increasing the former also leads to the increase of the latter. 



In expression (44), the term 6/a,3,f + 1) from (40) was neglected in comparison with 

6/a, 3, f ). This is a valid approximation for small \a\ and < f < 1; the reason why 

we are only interested in this regime is explained below. 

For < M < 1, 6/a,3,M)| as a function of the complex parameter a has a sharp 

peak at a = 0. Therefore, we may expect that the maximum of x (and r/) is obtained at the 

2 

minimum of |a| = ^ + ^ • For a lossless magnetic metamaterial slab, \a\ can be precisely 
zero when the perfect lens condition 

= -Q;/i„ (45) 

is met. In the point-dipole approximation, a lossless magnetic metamaterial slab can thus 
provide an infinite mutual inductance, leading to the perfect power transfer efficiency, rj = 1. 

In a real system, rj is limited by losses; its maximum is achieved by minimizing |a|, i.e. 
by requiring 

Re fix = —Re{a)fiy. (46) 

In this regime, a closed-form analytic expression for the mutual inductance can be obtained in 
the limit of small losses. Suppose that the magnetic loss tangents (Tx^z, introduced according 



to fix,z = ^^x.zi^ + j^x,z)i are both small, and /i^ = —fi^ v^yU^Z/i^. Note that we assume that 
both and //^ are negative. Then a ^ cr^/fJ^l, where a = {a^ + CTz)/'^- The other parameters 



in the expression (44) become approximately (to the 0(cr) accuracy): 

b = A/fil (47) 
h = A/fil (48) 
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To the same accuracy, the expression (44) becomes 

R2/Z2 X 



V 



R2/Z2 + ^^Li-b/a,3,w)l + X 

^<^l{-4/a^3,{l+v + w)/2) 

X = 7 w ^' (50) 

( t + ^<^L(-4/a^ 3, v)) (f + iM-^/a^ 3, w)) 

where we have introduced two parameters with the dimensions of impedance: 

^^'^ ^ ^°A(2^- ^^^^ 

Here, A = lixcjuj is the vacuum wavelength and = /igc is the free-space impedance. 

Before we find the optimum couphng regime, we can make two interesting observations 
about the behavior of Lerch function at large negative values of its first argument. First, 
one can show that 



lim -$i(-6/a, 3, u) 



00, < u < 1, , ^ 

(52) 



^ < 00, M > 1. 



Therefore, with arbitrarily small loss tangents, arbitrarily large enhancement of the mutual 
inductance L12 can be obtained when < m < 1, i.e. di + d2 < aD. 
Second, the slab contribution to the coil self-impedance 

L«oc-<l>^(-4/a^3,^;) (53) 
a 

diverges in the limit a — ?► whenever v < 1/2, i.e. when di < aD/2, but converges to a 
finite limit for any f > 1/2, in agreement with the findings of Refs.'^^'^. This means that the 
dipole source can be completely quenched by a lossless negative-/! slab in the superlensing 



regime (46) whenever the source is closer to the slab than the quenching distance, dg = aD/2. 
However, one should realize that with a finite loss, Ln remains finite for any di > 0; total 
quenching is possible only with a non-physical choice a = 0. 

The function ri{v,w) has a well-defined maximum, defining the optimum positions of 
the coils relative to the slab, as seen in Figure |2j The optimum positions depend on the 
resistances -Ri,2; a larger resistance Ri {R2) requires coil 1 (coil 2) to be closer to the slab's 
left (right) interface in order to achieve maximum performance. 

The maximum value of the function r], maximized over the distance variables {v, w), is a 
monotonically decreasing function of the loads Ri and i?2- It can be computed numerically 
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(a) (b) (c) 

FIG. 2. (Color online) Power transfer efficiency r/ as a function of the normalized distances v 
and w. (a) Ri/Zi = 1, R2/Z2 = 10, cr = 0.1. (b) Same as (a) but a = 0.01. (c) Ri/Zi = 10, 
R2/Z2 = 10, cj = 0.01. 
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FIG. 3. (Color online) Left: Power transfer efficiency r], maximized over normalized distances 
{VjW), as a function of the normalized resistance R2/Z2 for Ri/Zi = 1. For comparison, power 



transfer efficiency in free space - r]^""^ from equation (54 - is also shown as the black dotted curve. 
Right: The optimum normalized distances v (solid curves) and w (dashed curves) in the same 
configuration. Curve coloring corresponds to loss tangents a = 0.001 (blue), 0.01 (green), and 0.1 
(red). 

for any given loss tangent a. Figure [Sj^a) shows the maximum power efficiency for several 
choices of a. The optimum positions of the coils, fmax and Wmax, are shown in Figure |3](b). 
To find the regime in which a metamaterial slab offers significant improvement for the 



power transfer system, compare our result ( 50 ) with the no-slab case, for which 



1] = 7] 



X = X 



1 + x""" 

^vac ^vac 
I 1 2 

Ri R2 



(54) 
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where Z^f = Zi 2\aD=d- The expression (54) is small whenever either or both of the circuits 
are loaded with resistances exceeding the reference impedances Z^'^^- example, if both 
-Ri,2/^iT ratios are of order 10, x^"'^ is of order 10~^. For the same parameters, however. 



X given by expression (50) is of order 0.1 — 0.3 even with realistic loss tangents in the 
0.1 — 0.01 range. In this regime, the slab offers at least one order of magnitude improvement 
of the power transfer efficiency relative to the no-slab case. For even smaller loss tangents, 
efficiency r] from expression (50 ) can be close to 100%, as seen from Figure[3|^a). As expected, 
rj decreases monotonically with increasing loss tangent; it diverges slowly (logarithmically) 
in the limit cr — )■ 0. 



VI. THE EFFECT OF DIPOLE ORIENTATION 

So far, we have considered one orientation of the dipoles, namely, two dipoles parallel to 
the slab and to each other. Our formalism lets us consider two additional configurations: (a) 
two dipoles polarized in the ^-direction, and (b) one x-oriented and one z-oriented dipole. 
To consider these cases, our previous derivations must be complemented with a calculation 
of magnetic field of the 2;-oriented dipole. Following CheW^, the field of a magnetic dipole 
with moment rhi = rrtiz pointing in the positive 2;-direction is written as 

f oo jji 



H^ip, 0, ^) = ^ / dk,-^H^^\k,p)TTEe^'^^'-''\ (55) 

im k^ 

H,{p, cp,z) = ^-^ dkp-^H^^\k,p){l + RTEe-'^^')e^^^'\ (56) 

J-oo '^z 



in the transmitted {z > D) and refiected {—di < z < 0) regions, respectively. For the 
magnetic field at the locations of the two dipoles we obtain 

im k"^ 

H'r^\p = 0,z = D + d2) = ^--^ dkp-^TTEe'"^''', (57) 

4vr Jq kz 

Hr^ip = 0,z = -di) = ^ dk,'^RTEe'^''^'K (58) 

47r Jq kz 



These fields are different from the expressions (|29||30|) for the 

i^x ~ f^x case by a factor of 
two, but are otherwise the same. This means that L^^^ and L21 = L12 for the — 
case are twice the corresponding inductances in the — mx case. The coupling efficiency 



formula (50) suggests that x for the — case with load resistances Ri^2 is the same as x 



for the rux — rux case with resistances R'12 = -^1,2/2. Since the coupling efficiency increases 
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with decreasing loads, the — polarization provides a somewhat better coupling, with 
all other parameters fixed. The exact amount of additional enhancement can be determined 
from the diagram shown in Figure [Sj^a). 

Finally, we note that in our configuration with the two dipoles aligned to the z-axis, 
there is no coupling between rrix and rriz components: the axial magnetic field Hz of the 
x-polarized dipoles vanishes on the 2;-axis. 



VII. ON THE FEASIBILITY OF EFFICIENT COUPLING OF MAGNETIC 
DIPOLES THROUGH A DIELECTRIC-ONLY SLAB 

So far, we have been assuming that the coupling between the magnetic dipoles is mediated 
by the TE waves, which can be greatly enhanced by a negative-// slab; the contribution of 
the TM wave was neglected. Considering that metamaterials with electric-only response 
may offer lower loss tangents, it is interesting to see if a negative-e slab can compete with a 
negative-// slab. 

Consider a situation where only the TM waves are enhanced, such that the TM terms in 



expressions (29) and (30) dominate. This regime is achieved by choosing a dielectric-only 
slab with unit relative permeability and anisotropic relative permittivity tx^z < 0, and then 
tuning the real parts of e^^z to the superlensing condition, Ree^ = —Re{aTM)^v The mutual 



inductance (32) in this regime becomes 

/i0^1^2 



^21 



Air 



1 Ak^ 4 



(59) 



where = {(ye,x + ^e^z)!'^ is the effective dielectric loss tangent. The first term in (59) 



represents TE-wave mediated magnetic field transmission in free space. The second term in 



(59) dominates over the former whenever 



|*.(-;i,l,.™)»(^) . (60) 

Mathematically, it is possible to find a sufficiently small Oe for any < utm < 1 to satisfy 
this inequality, even with arbitrarily large A/rf ratio. The slab contribution to self-inductance 
is then dominated by the TM waves: 

(1) /ipA^fcg -2j 4 
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In this regime, power transfer efficiency looks formally the same as the equation (50), 
in which all quantities are replaced with their TM analogs, except that the characteristic 
impedances Zi 2 are now equal to 

Since the characteristic impedances set the scale for attainable load resistances -Ri,2, 
the dielectric slab can support only much smaller loads in the regime with reasonably high 
power transfer efficiency, rj > 0.1, in comparison with the magnetic slab. To achieve the 
same coupling efficiency with the same loads Ri^2 that the negative-/! slab can sustain, the 
loss tangent must be decreased to an exponentially small value, which is not feasible in 
practice. This finding is not applicable to two-dimensional (infinitely long) systems, where 
efficient coupling of two-dimensional point magnetic dipoles mediated by TM waves can be 
achieved with negative-e slab^. 



VIII. CONCLUSIONS 

The rigorous analysis presented in sections [Tl]jV] enables us to make the following conclu- 
sions about the metamaterial-based power relay system. 

First, for given frequency u, metamaterial slab thickness D, magnetic anisotropy ratio 
ote, and source and load resistances -Ri,2, there exists a global maximum of the power 
transfer efficiency 77 as a function of metamaterial permeability fiz and the positions of the 
coils di 2 relative to the slab. It is also possible to show that a global optimum exists with 
respect to variable parameters /i^, aTE and D when u, R12 and the total distance between 
the coils d = di + d2 + D are fixed. 

The two parameters - slab thickness D and anisotropy ratio ote ^ can leverage each 
Q^]^gjl30i35l Electromagnetically, the TE wave properties of an anisotropic slab of thickness 
D and anisotropy ote are indistinguishable from those of an isotropic slab of thickness 
D' = Dute- This statement can be derived using the invariance of Maxwell's equations 
with respect to coordinate transformations, a concept dubbed transformation optics in the 
recent year^. In practice, this means that the thickness of a metamaterial slab can be 
reduced without affecting performance, assuming that the magnetic loss tangent a can be 
kept constant while increasing the magnitude of the negative real part of 
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As a function of the useful load (-R2), optimum power transfer efficiency decreases slowly, 
much slower than the power law I/-R2 pertaining to the no-slab (vacuum) case. Consequently, 
there exists a threshold resistance Rt such that the metamaterial system loaded with R2 > 
Rt performs better than no-slab system. This threshold resistance Rt decreases with the 
metamaterial loss tangent a: less lossy metamaterials "beat" free space performance at 
smaller loads. 

Similar statements can be made about the efficiency as a function of the distance between 
the dipoles, d. Assuming that the anisotropy ratio is fixed, optimization of rj gives certain 
fixed ratios of the distances di 2 to the slab thickness D, and thus also the optimized ratio 
D/d. Since the power transfer efficiency depends on D (and, after optimization, indirectly 
on d) only through the ratios i?i,2/^i,2, where Zi^2 oc l/D^ oc 1/d^, Figure [i] can be used 
to determine how much rj decreases with increasing D. In the no-slab case, rj falls off with 
distance very quickly, as 1/d^. Again, the advantage of the metamaterial slab here is that 
it extends the range of distances d within which acceptable efficiency can be obtained. 

As typical for metamaterial superlenses, the performance can be theoretically arbitrarily 
good - in our case, the power transfer efficiency can be as close to 100% as needed - but 
only if arbitrarily small loss tangents can be implemented. The growth of the figure of merit 
with decreasing loss is very slow (typically, logarithmic). Nevertheless, even with practically 
attainable loss tangents a ~ 0.1 the metamaterial relay system can over-perform free space 
coupling efficiency by an order of magnitude or more, depending on the load resistance. 
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APPENDIX: CALCULATION OF REFLECTION AND TRANSMISSION CO- 
EFFICIENTS FOR A SLAB 

Consider a slab of homogeneous, uniaxially anisotropic medium occupying the space 
between < z < D, and a point magnetic dipole on the 2;-axis at z = —di, as shown in 
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FigjTj The spectral components of electric and magnetic fields can be derived from the 
and Ez fields of the TE and TM waves, whose evolution along the 2;-axis can be described 
as follows: 



-mi .,.2^(1), 



57r 



cos (j)k H{ (kpp) X < 



-di < z <0, 
0< z < D, 
z> D, 



(63) 



X02 



In the above expressions, 



XtE = -^J k'^Px/ l^z — A^xCi^O; 

^"p^x/^^z (^xf^xkf)- 



Xtm 



-di < 2: < 0, 
Q<z<D, (64) 
z> D. 



(65) 
(66) 
(67) 



From the continuity of Ey 



and Hx 

ox 



1 dEy 



juifix dz 



1 d'^H, 
fcp dxdz ' 



it follows that 



and are continuous. From these requirements and equations (|63|), one obtains four 



equations for the four unknowns Rte.Tte, ii.'Ti: 



e->^°'^>„(l + i?TE) =/ix(ti + ri). 



e-xodi^^(_l + /^^^) = ;^^^(-t, + n), 
XTi?(-tie-'^^^^ + rie^^^^) = -XoTte- 
The relevant solutions are 



(68) 
(69) 
(70) 
(71) 



- i) [exp(XTi?I^) - eM-XTED)] 

4?£exp(-Xo^i) 

+ £)'exp(xTi.i^) - - £)'exp(-xri./^) 
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Similarly, for TM polarization we use continuity of e^-Ez and with equations (64) to 



dz 

yield 



e-xo<ii,^(l + Rtm) = e.(pi + gi), (73) 

ex(pie-^-^^^ + gie>^-^^^) = e,T™, (74) 

e-^«"^Xo(-l + i?TA/) = Xtm{-Pi + gi), (75) 

XrM(-Pie-^™^ + gie^^*^^) = -XoTtm- (76) 



The coefficients of interest are 



- f ) [exp(xTM/^) - e^vi-XTAiD)] 
+ ^)'exp(xTAf/^) - - ^)%xp(-XTM/^) 



4XTMX0 (_ ^ ) 

^TM - Ta z r2 • {'') 

+ eMXTMD) - - exp(-XTM/^) 
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